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Abstract 

We introduce a method for classical trajectory calculations to simulate collisions between atoms 
and large rigid asymmetric-top molecules. Using this method, we investigate the formation of 
molecule-helium complexes in buffer-gas cooling experiments at a temperature of 6.5 K for molecules 
as large as naphthalene. Our calculations show that the mean lifetime of the naphthalene-helium 
quasi-bound collision complex is not long enough for the formation of stable clusters under the 
experimental conditions. Our results suggest that it may be possible to improve the efficiency of 
the production of cold molecules in buffer-gas cooling experiments by increasing the density of 
helium. In addition, we find that the shape of molecules is important for the collision dynamics 
when the vibrational motion of molecules is frozen. For some molecules, it is even more crucial 
than the number of accessible degrees of freedom. This indicates that by selecting molecules with 
suitable shape for buffer-gas cooling, it may be possible to cool molecules with a very large number 
of degrees of freedom. 
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I. INTRODUCTION 



Recent experimental work with ultracold molecules has generated a revolution in molec- 
ular physics [IHZ]- However, most of the experimental and theoretical studies in this field 
to date have focussed on diatomic or few-atom molecules [7HS]. An emerging direction in 
this field aims to extend the cooling techniques and study of ultracold diatomic molecules 
to large polyatomic molecules [10, llj. Cooling complex polyatomic molecules to cold and 
ultracold temperatures, if possible, will open up exciting opportunities for studying novel 
phenomena in physics and chemistry. For example, Kiipper et al have recently demonstrated 
how to manipulate the motion of polyatomic molecules in molecular beams using external 
fields jTU] . This may lead to the development of new experiments on precision measurements 
of large molecule collisions. One can also use intense laser fields [12j or static electric fields 
|13j to align cold molecules along a certain axis in the laboratory frame, which may provide 
a diversity of applications in molecular optics [2HT7], controlled chemistry [18], and spec- 
troscopic analysis 120] . Cooling ensembles of polyatomic molecules is also necessary for 
the research of cold chemistry [TJ. 

Among all cooling techniques, buffer-gas cooling is the most powerful and versatile 
method that produces molecules in translationally and internally cold states |2Tti23| . This 
technique relies on energy thermalization of molecules in a buffer gas (e.g., helium) main- 
tained at a cold temperature. Until recently, it was expected that polyatomic molecules, 
when placed in a cold buffer gas, will collect buffer gas atoms into clusters, which should make 
cooling polyatomic molecules difficult. However, Patterson et al have recently measured the 
thermalization rates of naphthalene (N) molecules in a cold helium gas [11] . Contrary to 
the expectation, they observed no clustering. They proposed an explanation based on the 
Lindemann mechanism [11] , Their model suggests that if the energy of the first excited 
vibrational state of the molecule is higher than the molecule-helium binding energy, there 
should be no clustering observed. However, this model excludes the role of rotational de- 
grees of freedom in cold collisions. Could the rotation induce a collision complex and thus 
trapping-clustering via a Lindemann mechanism? 

Here we undertake a detailed study of the collision dynamics between helium and two 
kinds of molecules, naphthalene and benzene, under conditions close to buffer-gas cooling 
experiments. At temperatures below 10 K, these molecules reside in the ground vibrational 
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state. We therefore treat the molecules as rigid. The rigid approximation can be justified 
semiclassically by noting that the vibrational modes in question are really ground state 
wavepackets represented by a mean position at rest. The conventional classical trajectory 
calculations for high energy dynamics allow all atoms in a molecule to move. [24j. These 
calculations, however, prove to be difficult to preserve the rigidity of molecules and are hard 
to implement for low temperature collisions excluding vibrational transitions. In addition, 
these methods meet with the zero point energy problem that may lead to unphysical energy 
localization in particular vibrational modes. Here, we adapt a molecular dynamics theory 
used in the simulation of polyatomic fluids to study low temperature collisions [2"5H2"T] . In 
our approach, the number of differential equations to solve is independent of the number of 
atoms in the molecule. Only the center-of-mass of molecules participates in the equations of 
the relative motion in a space-fixed (SF) frame. The rotational motion of the molecules is 
described separately by quaternion parameters in a body-fixed (BF) frame. We maintain the 
shape of the molecules by converting the coordinates of atoms in the molecules from the BF 
frame to the SF frame by a rotation matrix constructed using the quaternion parameters. 
The rigidity of the molecules is therefore preserved without introducing internal coordinates 
to constrain the bond length and angles of the molecules. 

Low energy collisions have been an interesting research topic for a long time |2"%l |2"U] . 
However, classical trajectory studies of collisions involving rigid molecules are restricted 
to small and symmetric-top molecules. This paper provides the first classical study of 
collision dynamics between a large rigid asymmetric-top molecule and an atom. The method 
introduced here can also find many other applications. 

II. METHODOLOGY 

We carry out classical trajectory calculations for naphthalene-helium and benzene-helium 
collisions at a temperature T = 6.5 K. In this section, we illustrate the methodology by 
naphthalene molecules. Naphthalene is a complex asymmetric-top molecule which has a 
large number of degrees of freedom. As shown in Figure. [TJ there are over one hundred 
non-degenerate rotational energy levels below 7 cm -1 . The complexity of the naphthalene 
molecules precludes the fully quantum calculation of the He-naphthalene collision dynamics. 
At low temperatures relevant for the buffer-gas cooling experiments, however, particles may 

3 



behave classically. The thermal de Broglie wavelength of helium at the temperature of 6.5 K 
is about 3.56 A. The mean interatomic distance in the buffer gas is around 136 A. This 
indicates that quantum effects in the collision dynamics are insignificant [30J. Moreover, 
a typical classical action associated to energy transfer between helium and naphthalene 
at this temperature is on the order of 10 cm -1 . It is much larger than the rotational 
energy splitting of naphthalene, as shown in Fig. [TJ suggesting that the dynamics under 
experimental conditions falls into the classical regime. We therefore employed a classical 
dynamics approach. In order to calibrate the accuracy of our calculations, we calculate the 
total thermal cross section for N— He collisions at T = 6.5 K. Our results are consistent with 
the experimental data. 

A. The system 

Fig. [2] demonstrates the coordinate system for collisions between naphthalene and helium 
at cold temperatures. The relative motion between naphthalene and helium is described 
in a space-fixed frame (x, y, z) with the origin (O) located at the center of mass of the 
collision complex N-He. The inter-particle distance R is the position vector of He relative 
to the center of mass of molecule N. The rotational motion of N is described using a body- 
fixed frame (x',y',z f ) with the origin positioned at the center of mass of N. The x',y', 
and z' axes correspond to the principal axes of the molecule. The rotational constants of 
naphthalene molecules used in this calculation are: A = 0.10405 cm -1 ; B = 0.04113 cm -1 ; 
C = 0.02948 cm -1 [3 lj . The orientation of the molecule N with respect to the space-fixed 
system is described by Euler angles. When we integrate the equations of rotational motion 
of molecules, we transform Euler angles into quaternion parameters [25-27J in order to 
eliminate the singularities in the equations of motion for the molecule. The approach is 
described in detail in the next section. The interaction potential surfaces between N and 
He are constructed based on the atom-bond pairwise additive representation developed by 
Pirani et al 132, 331. 
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B. Collision dynamics 



The relative motion of naphthalene and helium is described in a space-fixed Cartesian 
coordinate system. The motion of particles is determined by numerically integrating a 
system of classical equations of motion: 



dR p 

dt 



dp 

~dl 



r>SF 



(1) 
(2) 



where t and p are the time and the momentum for the N— He relative motion, respectively, 
and [i is the reduced mass. F SF is the interaction force between N and He, which is given 
by the first derivative of the interaction potentials V(R), 
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The orientation of naphthalene is described by quaternion parameters in a body-fixed 
frame, which are defined as follows: 
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where 9, <p, and if) are Euler angles. The rotational motion of naphthalene molecules is 
determined by a set of second order differential equations |27| : 



(8) 



where m = 0, 1, 2, 3; q m and q m are the first and the second time derivatives of q m , respec- 
tively; oj X) y^ z are the components of the angular velocity around the principal axes in the 
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body- fixed frame, and cu x ,y,z are their first derivatives with respect to t. In Equation [8j co x ,y,z 
can be obtained from the Euler's equations: 

1 
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where are the principal moments of inertia of the molecule and N xy z are the compo- 
nents of the torque in the body- fixed frame. We can calculate N xy z by transforming the 
torque in the space-fixed frame N X F z with a rotational matrix A |25l - l27] 



N BF = ^ N SF ; (12) 

defined as follows 



A 
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The torque is exerted on the center-of-mass of the molecule due to the interaction between 
the molecule and the helium atom. N is then given by: 

N SF = Rx (— F SF ). (14) 

In this calculation, we use the fourth-order Runge-Kutta method to solve coupled dif- 
ferential equations for both translational and rotational motions of the particles. The time 
stepsize for the integration is 5 x 1CT 3 ps. For low collision energies, e.g., < 1 cm -1 , the 
stepsize is chosen to be 1CT 5 ps. The total energy during collisions can be conserved up to 
three digits and the structure of naphthalene remains rigid throughout the simulation. The 
trajectories are terminated at large inter-particle separations (-Rfi na i) where the interaction 
potentials can be ignored. 

In this paper, we calculate two physical quantities: (i) the total cross sections for collisions 
between naphthalene and helium; (ii) the lifetime of the quasi-bound collision complex 
N— He*. We choose -Rfi na i to be 50 A in the calculation of (i) and 20 A for (ii). The 
lifetime of the collision complex N— He* is determined as the time between the first and 
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the last turning points. The rotational energy of naphthalene can be calculated from the 
final principal rotational velocities of naphthalene oo' xyz , which are given by the following 
equation 
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The trajectories for numerical calculations are generated with the initial conditions chosen 
as follows. The initial rotational energy of the molecule E TOt is determined by diagonalizing 
the matrix constructed with the hamiltonian for the rotational motion of a free asymmetric- 
top rigid rotor. We randomly sample the Euler angles and the magnitudes of u x , y ^ z with the 
constraint I x &x + ^y^y + = E™t to ensure that the trajectories cover the entire phase 
space jH]. The impact parameter hi is chosen quantum mechanically |34j : 

h = ^/I(r+T)h/P SF , (16) 

where / is the quantum number for the angular momentum describing the relative rotation 
between naphthalene and helium. For the calculation of (i), the value of the largest impact 
parameter 6™ ax is determined so that no trajectories have rotational energy transfer larger 
than the smallest energy gap between the rotational energy levels of naphthalene. For (ii), 
6™ ax is chosen to be large enough so that the lifetime of N— He* obtained from the calculations 
of the trajectories with hi > 6™ ax is zero. We carry out calculations of 4 x 10 5 trajectories for 
each collisional energy and the results were converged up to two digits. The total number 
of trajectories in this calculation is on the order of 10 7 . 

III. RESULTS 

In classical trajectory studies of collision dynamics of molecules, due to the accumulation 
of numerical errors, the conservation of the total energy of the collision complex and the 
rigidity of the molecules are generally difficult to maintain. In order to check the validity 
of our calculations, we calculate the total energy of the N-He collision complex and the 
C— C and C— H bond lengths of naphthalene as functions of time, as shown in Fig. [3] and 
Fig. [4], respectively. Despite the fluctuations of the total energy in the regime of short-range 

7 



interactions between naphthalene and helium, the total energy of the collision complex is 
intrinsically conserved to within three digits. The structure of molecules remains rigid 
throughout the simulation. 

At low temperatures, collisions between inert gas atoms and a complex molecule may 
induce energy transfer between translational and internal degrees of freedom. As a result, 
the atom may be trapped and move around the molecule for a long time leading to the 
formation of a long-lived quasi-bound collision complex. We observe some trajectories in our 
calculations where the exchange between relative translational energy and rotational energy 
of the molecule is very active, as demonstrated in Fig. [5j In this trajectory, the translational 
energy (black curve) oscillates significantly due to the interaction between naphthalene and 
helium and there are many rotational excitations and de-excitations of naphthalene (red 
dotted curve) during the collision process. At the end of the trajectory, the rotational energy 
of naphthalene is elevated from 4.75 cm -1 to 8.85 cm -1 ; correspondingly, the translational 
energy decreases from 6.00 cm -1 to 1.89 cm -1 . Fig.|6]shows for the same trajectory the inter- 
particle distance between naphthalene and helium as a function of time. The helium atom 
is trapped near the collision center from 24.00 ps to 74.25 ps where a long-lived collisional 
complex is formed. The lifetime of the complex in this trajectory is about 50.49 ps, which 
is determined as the time the particles spend from the first to the last turning points, as 
shown in Fig. [5] Fig. [7] displays for the same trajectory the z coordinates of the relative 
motion between helium and naphthaline as a function of its x and y coordinates. This graph 
depicts the trajectory in three dimensions. It demonstrates an orbital feature of the motion 
of the particles. 

In order to explore the role of the molecular shape in determining cold collisions, we 
calculate the average lifetime of the quasi-bound collision complexes B— He* and N— He* as 
a function of the collision energies, as shown in Fig. [8] by red curves and blue curves, respec- 
tively. Each point in this graph is the lifetime averaged over all trajectories with rotational 
energies of molecules below 20 cm -1 according to the Boltzmann distribution. Both curves 
display a decay of the lifetime as the collision energy increases. However, contrary to expec- 
tation, the lifetime of the B-He complex is longer than that of N-He even though benzene is 
lighter than naphthalene. Naphthalene is an asymmetric top while benzene is a symmetric 
top. We have established that the symmetry of benzene plays a twofold role in increasing 
the lifetime of the collision complex. First, it enhances the probability of the occurrence of 
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long trajectories. Based on our observations, long trajectories mainly occur when helium 
approaches the molecules near their edges. For a certain impact parameter which may lead 
to a long-lived collision complex, benzene allows for more approach directions than naphtha- 
lene that generate long trajectories. As a result, benzene molecules have a larger fraction of 
long trajectories than naphthalene. This mechanism is clearly demonstrated in Fig.[9j which 
shows the fraction of classical trajectories at a particular lifetime interval as a function of the 
lifetime of the B— He* (red) and N— He* (yellow) complexes at a collision energy of 2 cm -1 
for J = 0. Naphthalene has a bigger fraction of short trajectories (lifetime < 10 ps) than 
benzene, while benzene has a larger fraction of long trajectories. This difference becomes 
more significant as the lifetime of the collision complex increases. Second, molecular sym- 
metry may enhance the number of rotational periods in a molecule-helium complex, which 
leads to trapping of helium for a very long time. To illustrate this, we compare the lifetime 
of the longest trajectory for benzene and naphthalene. The longest lifetime for B— He* in 
our calculation is 1543 ps, while for N— He*, it is only 592 ps. These results suggest that the 
geometry of molecules plays a significant role in determining the dynamics of cold collisions. 

The long-lived quasi-bound collision complex, as demonstrated in Fig. |5j |6j and [7| could 
result in a stable molecule-He complex if the lifetime of the complex is long enough to allow 
subsequent collisions with helium. Here, we average the lifetime of the N— He* and B— He* 
complexes over the collision energy according to the Maxwell-Boltzmann distribution and 
compare it with the mean time a helium atom needs to collide with a naphthalene molecule 
7~He- The mean lifetime of the N— He* and B— He* complexes is 7.9 ps and 9.9 ps, respectively. 
THe can be estimated as an inverse of the number of collisions M one naphthalene molecule 
experiences in a helium gaseous bath per second, which is given by 



A^AWTN-tM/-^^- ( 17 ) 

^MN-He 



In Eq. 17, An e = 4 x 10 17 cm -3 is the density of helium atoms used in the experiments 
[11 j . o"N-He ~ 1CT 13 cm 2 is the total thermal cross section for collisions between naphthalene 
and helium |llj . k-Q is the Boltzmann factor, and T = 6.5 K is the temperature of the 



naphthalene-helium gaseous mixture. Using Eq. 17, we obtain M ~ 7.24 x 10 s s _1 and 
correspondingly, th g ~ 1380 ps. This time is two orders of magnitude longer than the 
calculated average lifetime of the complexes. The time a helium atom needs to collide with 
the quasi-bound collision complex should be even longer than Tjj e . This indicates that the 



mean lifetime of the quasi-bound collision complex N— He* is too short to allow subsequent 
collisions with a second He to form a stable molecular complex. As a consequence, no 
clustering occurs for naphthalene in a low-temperature He gas. This suggests that the density 
of the helium gaseous bath used in the buffer-gas cooling experiment could be increased up 
to 10 19 cm -3 , which may decrease the loss of the produced molecules during the cooling 
process and enhance the efficiency of cooling. 

We also calculate the thermal total cross section for collisions between naphthalene and 
helium. The cross section is 1.19 x 1CT 13 cm 2 , which is consistent with experimental obser- 
vation (Tn-Hc ~ 10~ 13 cm 2 [llj. 

IV. SUMMARY 

We have introduced a method for classical trajectory calculations of collisions between 
atoms and large rigid asymmetric-top molecules. We explore the collision dynamics between 
helium and large polyatomic molecules at low temperatures. In particular, we investigate 
the formation of naphthalene-helium and benzene-helium complexes in buffer-gas cooling 
experiments. Our finding provides an explanation for the experimental observation. Clus- 
tering would occur if helium atoms formed long-lived quasi-bound collision complexes with 
naphthalene molecules. Subsequent collisions with helium would result in the formation 
of stable naphthalene-helium complexes. Our calculations show that the mean lifetime of 
the collision complex N— He* is too short to allow subsequent collisions. As a consequence, 
no clustering occurs under the experimental condition. In the experiment, the density of 
helium atoms is about 4 x 10 17 cm -3 . Our calculations show that this density could be 
increased to 10 19 cm -3 for the buffer-gas cooling of naphthalene molecules. In buffer-gas 
cooling, the density of helium determines the time scale for the thermalization of transla- 
tional degrees of freedom. This time scale should be short enough to guarantee that the 
target molecules can be cooled to the desirable temperature before they are lost on the wall 
of the cell. Increasing the helium density therefore could potentially improve the efficiency 
of the production of cold molecules. In addition, we calculate the lifetime of the quasi-bound 
complex for collisions between benzene and helium at the temperature of 6.5 K. Contrary 
to the expectation, the mean lifetime of B— He* is longer than that of N— He*. This is due 
to the symmetric structure of benzene. Our results indicate that the shape of molecules is 

10 



important in determining the dynamics of cold collisions, in some cases, even more crucial 
than the number of accessible degrees of freedom. This suggests that by selecting molecules 
with suitable shape for buffer-gas cooling, one could cool molecules with a very large number 
of degrees of freedom. 
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FIG. 1: The rotational energy levels of naphthalene below 7 cm 1 . 
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FIG. 2: The coordinate system for collisions between naphthalene (N) and helium at cold tem- 
peratures. The relative motion of N and He is described in a space-fixed frame (x, y, z) with the 
origin (O) located at the center of mass of the collision complex N-He. The inter-particle distance 
R is the position vector of He relative to the center of mass of the molecule N. The motion of N is 
described using a body-fixed frame (V, y' , z') with the origin positioned in the center of mass of N. 
The x',y', and z' axes correspond to the principal axes of the molecule. 
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FIG. 3: The total energy of the N-He collision complex as a function of time at the collision energy 
of 6 cm" 1 . The energy of the system is intrinsically conserved up to within three digits. 
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FIG. 4: The C— C and C— H bond lengths of naphthalene as functions of time during the collision 
process. The rigidity of Naphthalene is fully preserved. 
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FIG. 5: The translational energy (black curve), rotational energy (red dotted curve), and interaction 
potential energy (blue curve) as functions of time during a collision process between naphthalene 
and helium. 
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FIG. 6: The inter-particle distance between naphthalene and helium as a function of time. The 
helium atom is trapped near the collision center from 24.00 ps to 74.25 ps. 
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FIG. 7: The z coordinates of the relative motion between helium and naphthalene as a function of 
its x and y coordinates for a long classical trajectory. 
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FIG. 8: The average lifetime of the quasi-bound collision complex of N— He* (red curve) and B 
(blue curve) as functions of collision energy. 
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FIG. 9: The fraction of classical trajectories at a particular lifetime interval as a function of the 
lifetime of the B— He* (red) and N— He* (yellow) complexes at a collision energy of 2 cm -1 for 
J = 0. Naphthalene has a bigger fraction of short trajectories (lifetime < 10 ps) than benzene, 
while benzene has a larger fraction of long trajectories. This difference becomes more significant as 
the lifetime of the collision complex increases. 
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